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Abstract 

We present results from the analytic calculation of ttH hadronic production at Next-to-Leading 
Order in QCD interfaced with parton-shower Monte Carlo event generators in the POWHEG BOX 
framework. We consider kinematic distributions of the top quark and Higgs boson at the 8 TeV 
Large Hadron Collider and study the theoretical uncertainties due to specific choices of renor¬ 
malization/factorization scales and parton-showering algorithms, namely PYTHIA and HERWIG. The 
importance of spin-correlations in the production and decay stages of a top/antitop quark is dis¬ 
cussed on the example of kinematic distributions of leptons originating from the top/antitop decays. 
The corresponding code is now part of the public release of the POWHEG BOX. 
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I. INTRODUCTION 


The associated production of a Higgs boson with a pair of top (t) and antitop quarks 
[i) has been receiving increased attention after the discovery of a Higgs boson (H) at the 
CERN Large Hadron Colllider (LHC) [U [2] . Although rare in nature, this production mode 
gives direct access to the Higgs-boson coupling to the heaviest elementary fermion, the top 
quark, and provides important constraints on extensions of the Standard Model (SM) when 
combined with the indirect determination of the same coupling in Higgs-boson production 
via loop-induced gluon-gluon fusion. Experimental studies of Higgs properties have already 
provided strong constraints in Run 1 [3HS], and will be one of the main focuses of Run 2 of 
the LHC. 

Recent experimental analyses have explored the possibility of detecting ttH production 
when the Higgs boson decays via H ^ bb [6H9], iL —)■ 77 [T01 - IT2] . and H —)■ ZZ*, WW* [T3] . 
or in a combination of all the previous modes |T3]. Apart from overcoming challenging 
background difficulties, it is important first of all to have the signal fully under control, 
including the decay of the hnal-state particles. With this respect, the theoretical systematic 
error has to be carefully investigated and reduced, if possible, by adding both QCD and 
electroweak higher-order corrections consistently interfaced with existing parton-shower (PS) 
Monte Carlo codes, e.g. PYTHIA [151 HH], HERWIG [HI HE] or SHERPA |I9], to match the 
complexity of real events. 

The parton-level O^a^) or next-to-leading order (NLO) QCD cross section for tiH produc¬ 
tion has been Erst calculated in Refs. jEOHES] and subsequently conhrmed by various collab¬ 
orations using one-loop computational tools such as aMCONLO [26l l27j and PowHel |2H1 129] . 
More recently the 0(asa'^) contribution to the parton level cross section has also been 
calculated including [30] or omitting [31] QED corrections. Effects of parton shower and 
hadronization on the NLO-QCD parton-level cross sections have been studied both in the 
aMCSNLO framework ^B] and in a private implementation of the POWHEG method by the 
PowHel collaboration [2H]. The two approaches have been compared and found in good 
agreement in a study performed within the LHC Higgs Cross Section Working Group [32] . 
More recently, the original NLO-QCD calculation of Refs. [221 - 125] has been interfaced with 
SHERPA [19] and has become part of the SHERPA-2.0.0 release. We have now interfaced the 
same routines with the POWHEG BOX [E3HSS] including the decays of the top/antitop quarks 
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and the Higgs boson. The top-qnark spin-correlations have been taken into acconnt in an 
approximate way, similar to what has been done in the aMCONLO [Sn] and PowHel [2H] imple¬ 
mentations, to allow for studies based on top-quark spin-polarization effects |3SHSB] • A hrst 
comparison of this interface with PowHel and SHERPA, that did not include spin-correlation 
effects, was presented in the context of the Les Houches 2013 Workshop [39]. After more 
detailed cross checks and further improvements, we are now making our implementation 
public through the POWHEG BOX website, http://powhegbox.mib.infn.it/, 

In view of the dedicated effort on tiH analyses in both the ATLAS and CMS collabora¬ 
tions, in this paper we would like to present the details of our implementation and take the 
opportunity to address some of the issues that have emerged in recent experimental studies. 
In particular, we will investigate the dependence of theoretical predictions on the choice of 
a static or dynamical renormalization/factorization scale, the possibly different behavior of 
the interface with PYTHIA/HERWIG, and the effect of top-quark spin-correlations on various 
kinematic distributions of the decay products of the ttH hnal state. Even though the back¬ 
ground to tiH represents a major hurdle which needs to be overcome in order to properly 
measure this very important production channel, a very accurate control of the signal is still 
the hrst necessary step, and the possibility for the experimental community to have access 
to several public tools is very valuable. Our implementation answers this need and provides 
the original NLO-QCD analytic calculation |22] - l25] in the same framework (POWHEG BOX in 
this case) as other processes that enter the tiH studies, from ti production [101 to single-top 
production [H] or Higgs production in gluon-gluon fusion [121 SO] , offering a fully consistent 
alternative to analogous studies in the MadGraph5_aMC@NL0 [B] framework. 

II. IMPLEMENTATION 

The implementation of a new process in the framework of the POWHEG BOX requires a list 
of all independent havor structures for the Born and the real-emission contributions, the 
Born and the real-emission amplitudes squared, the hnite parts of the virtual amplitudes 
interfered with the Born amplitude, the spin- and color correlated amplitudes squared, and 
a parametrization of the phase space for the Born process. The tree-level amplitudes can 
be generated with the help of the build tool based on MadGraph 4 [l5l Il6] that is available 
in the POWHEG BOX. For pp —)■ tiH, we use the virtual amplitudes of [221125] and adapt them 
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to the format required by the POWHEG BOX. In order to validate our implementation, we 
have checked that the results of [21112S] are fully reproduced by the NLO-QCD mode of the 
POWHEG BOX. 

The hxed-order calculation of [221125] provides NLO-QCD corrections to the production 
of an on-shell ttH state at a hadron collider. If implemented in a multi-purpose Monte-Carlo 
program like PYTHIA, decays of the scalar Higgs boson can easily be simulated. In principle, 
also the decays of the top quarks can be taken care of externally. However, merely combining 
the on-shell calculation for pp —)■ tiH with a separate, spin-averaged simulation of the top- 
quark decays by the shower Monte Carlo results in a loss of information on correlations 
between the production and the decay stages. 

A method for overcoming this limitation was proposed in ffzi and has by now been applied 
in several POWHEG BOX implementations of processes with unstable particles [m SHI SS]. 
The basic idea of this approach is to hrst produce on-shell tiH events including NLO-QCD 
corrections matched to a parton-shower generator (NLO-pPS) via the POWHEG framework, 
and subsequently simulate the decays according to a distribution determined by matrix 
elements for H{t ——)■ M~9) and H{t ^ bi~^iy){i —)■ bi~P) -P jet hnal states that 
retain correlations in the top-quark decays, both at leading order and in the real-emission 
contributions. In addition to an improved description of spin correlations between the 
production and decay stages, this method allows for a moderate reshuffling of momenta 
in such a way that the virtualities of the top quarks are distributed according to Breit- 
Wigner shapes. For the actual implementation of this procedure we follow closely what 
has been done for the related case of ti -P jet production ra. and refer the reader to that 
reference for further technical details. The code version we are providing in the POWHEG BOX 
allows the user to activate this feature for an improved description of top-quark decays via 
a switch in the input hie. Alternatively, the code can generate on-shell events for pp —tiH 
that are subsequently decayed externally, for instance via PYTHIA. 


III. RESULTS 

The results presented in this section have been obtained using a prototype setup where 
we consider the LHC operating at a/s = 8 TeV, with the top quark and Higgs boson masses 
chosen to be, respectively, rrit = 172.5 GeV, and mn = 125 GeV. After having validated 
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the POWHEG BOX implementation against the parton-level NLO-QCD results of Refs. 

1^ . we have focused on some of the issues that have been raised in recent experimental 
analyses concerning differences encountered in the behavior of several distributions when 
interfacing the NLO-QCD calculation of pp — ttH with either PYTHIA or HERWIG. With 
this respect, results are presented using the PYTHIA-6.4.25 [15], PYTHIA-8.183 [16], and 
HERWIG-6.5.10 [18] releases of the aforementioned parton-shower algorithms. We note that, 
in principle, POWHEG requires a transverse-momentum ordered parton-shower program, such 
as PYTHIA. In order to correct for unwanted radiation effects that occur in angular-ordered 
parton shower programs, such as HERWIG, a vetoed-truncated shower has to be included. 
Since the public version of HERWIG does not provide a truncated shower, we neglect these 
contributions and consider them as part of the systematic uncertainties of our results. In all 
our numerical results, underlying event and multi-parton interactions are disregarded. We 
note that more specihc studies should be developed in the context of dedicated experimental 
analyses and we facilitate them by making our codes public through the POWHEG BOX. 

In this section we choose to illustrate the reach of our implementations by showing a 
sample of distributions for the main hnal-state objects (t, i, and H) and their decay products, 
obtained by interfacing the NLO-QCD calculation with both PYTHIA and HERWIG. Indeed for 
our numerical analysis, we only consider the decay products of the top and antitop quarks. 
Decays of the scalar Higgs boson (no spin-correlation involved) can easily be simulated with 
the external parton-shower Monte Carlo, and we provide a respective option in the public 
version of our POWHEG BOX implementation. We study the dependence of distributions on 
the choice of either a hxed or a dynamical renormalization/factorization scale. We do not 
take up a study of other issues, like the dependence on intrinsic parton-shower scales, or the 
residual uncertainty due to the choice of Parton Distribution Functions (PDF). These studies 
are important and should be performed in the context of specihc experimental analyses. All 
our results are presented using the CTlOnlo [SO] set of PDF as implemented in the LHAPDF 
library lai with the number of light quark havors set to Nf = 5 and the strong coupling at 
NLO QCD determined by as{Mz) = 0.118. 

To explore the issue of scale dependence, we set the renormalization and factorization 
scales equal to each other and studied both the case of a hxed and a dynamical scale. When 
using a hxed scale we have chosen /iq = rrit + mu/‘2^ as central scale, while we have chosen 
fiQ = {mT{t)mT{i)mT{H)y^^ as central value when working with a dynamical scale (where 
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FIG. 1; The pT (left) and t] (right) distributions of the Higgs boson at NLO-QCD with no 
parton shower (solid, black), and with parton shower as obtained through PDWHEG+PYTHIA6 
(long-dashed, red), P0WHEG-I-PYTHIA8 (short-dashed, orange), and POWHEG-I-HERWIG (dot-dashed, 
blue) respectively, for a fixed-scale choice (see text). The lower panels show the ra¬ 
tios: R = cicr(NLO)/d(T(PYTHIA6) (black), R = ficr(HERWIG)/dcj(PYTHIA6) (red), and R = 
(icj(HERWIG)/(icr(PYTHIA8) (blue). The error bars indicate the statistical uncertainties of the Monte- 
Carlo integration. 


mT{P) denotes the transverse mass of particle P). The renormalization/factorization scale 
dependence of the results is assessed by arbitrarily varying such scales by factors of = 1/2 
and 2 about the central value po. As default choice, we opt for the hxed scale and use 
Rr = Rf = 'fTT't + ^h/‘^ unless otherwise specihed. 

To study the impact of different parton-shower generators on the NLO-QCD parton-level 
results, we present in Figs. [TJj^the comparison between the parton-level NLO-QCD results 
and the corresponding results obtained upon matching with either PYTHIA6, PYTHIAS, or 
HERWIG as implemented in the POWHEG BOX framework. In Figs. [Tjj^we illustrate the impact 
of the parton shower on the hxed-order NLO-QCD results for the transverse-momentum and 
rapidity distributions of the Higgs boson and the top quark, respectively. Parton-shower 
effects do not cause large changes compared to the NLO results of these distributions, and 
differences between the PYTHIA and the HERWIG implementations are small. We noticed 
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FIG. 2: The pT (left) and r] (right) distributions of the top quark at NLO-QCD with no 
parton shower (solid, black), and with parton shower as obtained through PDWHEG+PYTHIA6 
(long-dashed, red), P0WHEG-I-PYTHIA8 (short-dashed, orange), and POWHEG-I-HERWIG (dot-dashed, 
blue) respectively, for a fixed-scale choice (see text). The lower panels show the ra¬ 
tios: R = dcr(NLO)/d(T(PYTHIA6) (black), R = ficr(HERWIG)/dcj(PYTHIA6) (red), and R = 
(icj(HERWIG)/(icr(PYTHIA8) (blue). The error bars indicate the statistical uncertainties of the Monte- 
Carlo integration. 


however a systematic enhancement of the low pt region in HERWIG with respect to PYTHIA6, 
while PYTHIAS seems to have a better agreement with HERWIG over the entire pr spectrum. 
The right panel of Fig. [^shows the transverse momentum of the ttH system. This observable 
is entirely due to real-radiation contributions in the NLO-QCD corrections and parton- 
shower effects. At leading order, the transverse momentum of the ttH system is zero, 
because of momentum conservation. If real-emission contributions are taken into account, 
as in the hxed-order NLO-QCD calculation, this observable exhibits a divergence as the 
transverse momentum of the entire system approaches zero. In the NLO-pPS result, this 
behavior is tamed by a Sudakov factor. The POWHEG-(-PYTHIA and POWHEG-(-HERWIG results 
are thus much better behaved in the region of low pT{ttH) that correspond to the emission of 
a soft jet, and they are compatible over the full pt spectrum. We could present distributions 
for jet observables as well, starting from the px and p of the hardest or next-to-hardest jets. 
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FIG. 3; The pT distributions of the tt pair (left) and of the ttH system (right) at NLO-QCD with 
no parton shower (solid, black), and with parton shower as obtained through PDWHEG+PYTHIA6 
(long-dashed, red), P0WHEG-I-PYTHIA8 (short-dashed, orange), and POWHEG-I-HERWIG (dot-dashed, 
blue) respectively, for a fixed-scale choice (see text). The lower panels show the ra¬ 
tios: R = cfcr(NLO)/d(T(PYTHIA6) (black), R = (icr(HERWIG)/do-(PYTHIA6) (red), and R = 
(icj(HERWIG)/dcr(PYTHIA8) (blue). The error bars indicate the statistical uncertainties of the Monte- 
Carlo integration. 


but they would not help us to judge the impact of the parton shower if we did not set 
in place, at the same time, more specihc selection cuts aimed at distinguishing the hrst 
emission described by the hard matrix elements or the top-quark decays from the following 
emissions coming from the parton shower, as well as cuts and vetoes aimed at distinguishing 
light jets from heavy-flavor jets. Since these sort of requirements only make sense in the 
context of dedicated experimental analyses, we refrain from considering specihc jet variables 
and present only results for observables of systems, like the ttH system, that carry a clear 
imprint of the jet from the hrst QCD emission, and allow us to verify expected ehects like 
the one due to the Sudakov factor. 

In order to assess the theoretical uncertainties associated with the choice of renormaliza¬ 
tion and factorization scale, we have computed the previously considered distributions for 
diherent choices of scale as previously explained. In particular, Figs. |4]j^are obtained using 
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FIG. 4: The pT (left) and rj (right) distributions of the Higgs boson obtained with PDWHEG+PYTHIA6 
for fixed (fix) and dynamical (dyn) renormalization/factorization scales. The lower panels show 
the respective ratios R = da{^po)/da{iJ,o) for ^ = (0.5; 2). 


our P0WHEG+PYTHIA6 implementation. More specifically, Figs. and show the transverse 
momentum and rapidity distributions of the Higgs boson and the top quark, respectively. 
The scale dependence of the results is considerable, amounting to more than ±10% in some 
regions of phase space. Using a dynamical rather than a hxed scale helps in slightly reducing 
the scale uncertainty of the NLO±PS results in all observables. 

A similar behavior can be observed in the transverse momentum distributions of the 
ti pair and the tiH system, c.f. Fig. Since the latter observable vanishes at leading order, 
we expect it to be plagued by larger uncertainties than distributions that are genuinely 
described at NLO accuracy. 

We finally illustrate the impact of spin-correlation effects. As pointed out, for instance in 
Ref. [3H] , exploiting polarization effects in the top-quark decays can appreciably improve the 
sensitivity of the LHC in the tiH channel. At the same time, spin-correlation observables 
can probe the top-Higgs coupling in tiH production at the LHC, as studied for instance in 
Ref. [52] • To allow for the simulation of spin correlations in the decays of the top quarks our 
POWHEG BOX implementation resorts to the prescription of Ref. [17], as explained in Sec. 
NLO-QCD corrections are considered for the tiH production process only. Improving on 
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FIG. 5; The pT (left) and p (right) distributions of the top quark obtained with P0WHEG+PYTHIA6 
for fixed {fix) and dynamical {dyn) renormalization/factorization scales. The lower panels show 
the respective ratios R = da{f,po)/da{iJ,o) for ^ = (0.5; 2). 
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FIG. 6: The pT distributions of the tt pair (left) and of the ttH system (right) obtained with 
P0WHEG+PYTHIA6 for fixed {fix) and dynamic {dyn) renormalization/factorization scales. The 
lower panels show the respective ratios R = da{(,po)/da{po) for ^ = (0.5; 2). 
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FIG. 7: Normalized transverse-momentum distribution of the hardest lepton (left) and az¬ 
imuthal angle separation of the two hardest oppositely charged leptons (right) obtained with 
P0WHEG-I-PYTHIA6 with (blue lines) and without (red lines) spin correlations in the decays of the 
top quarks, for a dynamic-scale choice (see text). The lower panels show the ratios of the results 
without and with spin correlations. 


the accuracy of the decay process would require more advanced techniques, such as the 
method that has recently been presented for top quark pair production at the LHC [5B] . 
For processes with additional particles in the hnal state, such as ttH production, such a 
procedure has never been applied so far and we restrict ourselves to the more approximated 
description of the decay process for the time being. To explore the impact of spin-correlation 
effects on experimentally accessible observables, we consider the leptonic decay modes of the 
top quarks and focus on final states with at least two oppositely charged leptons with 

PT,e> 20 GeV, \ye\ < 2.5 , (1) 

that are well-separated by ARij > 0.4 from identihed jets with ptj > 20 GeV and \yj\ < 4.5. 
As illustrated by Fig. the transverse-momentum distribution of a hard lepton barely 
depends on the treatment of the top-quark decays. Predictions for angular correlations 
between the two hardest oppositely-charged leptons, however, exhibit a sizable modification 
when the simulation of the top-quark decays is improved. 
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IV. CONCLUSIONS 


In this work we have presented a new, pnblic implementation of tiH production at the 
LHC at NLO-QCD matched to parton showers in the framework of the POWHEG BOX. In order 
to illustrate the reach of the implementation, we discussed several sources of theoretical 
uncertainties, in particular scale dependencies and differences due to the specific shower 
Monte Carlo program used. We found that, within the limits of such a general study, for 
most observables differences between predictions obtained with PYTHIA6, with PYTHIAS, and 
with HERWIG are small. We remind the reader that the public version of HERWIG that we 
have used does not have any vetoed truncated showers and therefore the comparison between 
PYTHIA and HERWIG will have to be further investigated in particular in analyses that include 
jet observables. Scale uncertainties are generally non-negligible and can amount to more than 
±10% in some regions of phase space. The uncertainties obtained using a dynamical and 
a hxed scale are compatible, although typically they are smaller when using a dynamical 
scale. For observables involving the decay products of the top quark, we recommend using 
a code version that does take into account spin correlations between the production and 
decay stages. 
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